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1.0  INTRODUCTION 


The  purpose  of  this  report  is  to  describe  the  results  of  studies  to 
investigate  the  application  of  standard  laser  velocimeter  (LV)  techniques 
for  measurements  in  turbulent  flows.  Turbulence  is  characterized  by 
the  excursions  of  the  instantaneous  velocity  at  a  point  about  a  mean  flow 
velocity.  Both  the  frequency  and  the  magnitude  of  these  excursions  are 
of  interest  to  the  aerodynamicist.  Studies  are  under  way  by  several 
organizations  (Refs.  1,  2,  and  3)  involving  LV  measurements  of  the 
power  spectral  density,  shear  stress,  and  other  parameters.  The 
scope  of  this  project  has  been  limited  to  measuring  the  magnitude  of  the 
turbulence,  the  turbulent  intensity. 

The  traditional  aerodynamic  diagnostic  tools  such  as  the  pitot  probe 
are  well  suited  with  their  time -averaging  aspects  for  the  measurement 
of  the  average  or  mean  velocity,  but  only  the  hot-wire  anemometer  has 
been  successful  in  measuring  turbulence  parameters.  The  hot  wire  pro¬ 
vides  an  indirect  measurement  of  turbulence  by  analyzing  the  heat- 
transfer  rate  of  the  wire  as  affected  by  convective  cooling.  It  is  extremely 
delicate,  requiring  extreme  care  in  use  and  tedious  calibration  procedures. 
The  LV  offers  many  advantages  over  the  hot  wire.  It  is  nonperturbing  to 
the  flow  of  interest,  and  it  is  unaffected  by  harsh  environments  such  as 
extreme  heat.  It  offers  very  good  spatial  resolution  with  typical  probe 
volumes  .as  small  as  2  x  10-5  cm^  and  allows  dose  proximity  to  models 
of  interest  for  boundary -layer  studies,  etc.  Its  outstanding  point  in  this 
study  is  the  inherent  similarity  of  the  LV  method  to  the  parameter  of 
interest.  The  LV  makes  many  individual  "instantaneous"  determinations 
of  velocity  from  which  a  mean  velocity  is  found.  Thus,  the  mean  velocity 
and  the  magnitude  of  the  excursions  about  the  mean  are  directly  mea¬ 
sured  by  the  LV. 

The  studies  conducted  under  this  program  were  intended  to  provide 
insight  into  known  problem  areas  concerning  LV  applications  and  to  un¬ 
cover  and  define  other  problems  through  the  reduction  of  typical  LV  data 
to  find  the  turbulent  intensity. '  Data  reduction  techniques  were  studied  in 
depth.  A  comparison  of  experimental  and  analytical  determinations  of 
accuracy  versus  sample  size  is  presented  for  mean  velocity  and  turbulent 
intensity.  It  is  hoped  that  direction  is  provided  for  implementation  of  the 
LV  technique  to  reliable  turbulence  measurements. 
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2.0  THE  LV  TECHNIQUE 


The  LV  technique  is  based  on  the  Doppler  frequency  shift  of  radia¬ 
tion  emitted  or  scattered  from  a  moving  source. 


The  shift  (or  difference)  in  frequency  between  radiation  incident  on 
a  moving  particle  and  the  radiation  received  after  being  reflected  or 
scattered  from  that  particle  is  expressed  as 


fr  u  V  /  A  A  \ 

s  -  f«  =  r  *  (es~ e* 

o 


(2.1) 


where 

fj  =  Doppler  frequency  shift 
fs  =  frequency  of  scattered  radiation 
fQ  =  frequency  of  incident  radiation 
n  =  index  of  refraction  of  the  medium 
"v  =  scatter  source  velocity  vector 
Aq  =  frce-space  wavelength  of  the  incident  radiation 

es  =  unit  vector  in  the  direction  of  scatter 

e;  =  unit  vector  in  the  direction  of  incident  radiation 


Equation  (2.  1),  derived  in  Ref.  4,  shows  that  the  measurable  fre¬ 
quency  difference  is  dependent  on  the  scatter  source  velocity  and  the 
directions  of  incident  and  scattered  radiation.  In  practice,  it  is  con¬ 
venient  to  use  a  particular  geometric  configuration  referred  to  as  a 
"dual  scatter"  system,  in  which  two  incident  beams  of  radiation  are 
simultaneously  scattered  from  the  same  particle  in  a  common  direction, 
es  (Fig.  1).  From  Eq.  (2. 1)  the  radiation  scattered  from  the  first  beam 
is 


fo  +  r  •  w.-V 

O  1 


(2.2) 


Similarly, 


t  ^ V  / A  /v  \ 

=  fo  +  -  *  <es-eJ 
o  * 


(2.3) 


The  frequency  difference,  measurable  by  the  optical  heterodyning 
technique,  is  then 


f 


(8;  -  ) 
'2  '1 


(2.4) 
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The  measured  frequency  difference  is  now  independent  of  the  direction 
of  observation  (es),  which  gives  great  versatility  to  this  arrangement. 
Further  convenience  is  added  in  the  form 


f 


d 


eA 


(2.5) 


where  9  is  the  angle  between  the  two  incident  beams  and  ex  is  the  unit 
vector  perpendicular  to  the  mean  illuminating  direction.  The  angle  9 
can  be  measured  quite  precisely,  thus  allowing  an  accurate  direct  mea¬ 
surement  of  scatter  source  velocity. 


Observer 


Figure  1.  Dual  scatter  technique. 

The  same  results  are  obtained,  but  more  insight  into  the  physical 
mechanism  is  acquired  by  looking  at  the  spatial  interference  pattern 
from  the  two  coherent  illuminating  beams  of  the  dual  scatter  system. 

A  set  of  bright  and  dark  fringes  is  set  up  where  the  two  incident  beams 
cross,  as  shown  in  Fig.  2.  The  input  waves  are  indicated  by  lines 
representing  the  high  intensity  crests  one  wavelength  apart.  Construc¬ 
tive  interference  of  the  two  produces  the  fringes,  bright  where  the  two 
wave  crests  are  superimposed,  and  dark  at  the  superposition  of  a 
crest  and  a  trough.  The  spacing  between  fringes  is  seen  to  be 


2  sin  d/2 
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A  particle  moving  across  the  fringe  pattern  will  scatter  light  with  in¬ 
tensity  varying  as  the  fringe  pattern  intensity  and  with  a  frequency 
directly  proportional  to  the  component  of  its  velocity  normal  to  the 
fringes.  The  period  of  this  a-c  signal  is  the  time  required  to  travel  a 
known  distance  (the  fringe  spacing) .  Fundamentally,  velocity  equals 
distance  divided  by  time,  or 

-*■  A  A  1 

v  •  e  =  -  x  — 

n  2  sin  0/2  T  (OR 


or  as  before,  since  the  signal  frequency  (f)  equals  1/r  and  the  wave¬ 
length  in  a  medium  of  index  of  refraction  n  is  A0/n,  then 


f  = 


sin  0/2 
— i -  nv 


(2.7) 


where  en  is  a  unit  vector  normal  to  the  fringes.  It  is  obvious  from 
Fig.  2  that  en  =  e±. 


Many  arrangements  using  this  basic  method  are  possible.  A  typical 
dual  scatter  system  is  schematically  shown  in  Fig.  3.  A  signal  burst  of 
velocity  information  is  received  each  time  a  scattering  source  passes 
through  the  fringe  region.  Eight  cycles  are  used  to  determine  the  period 
of  the  signal  by  the  Data  Processor  (DP)  developed  at  AEDC  (Refs.  4 
and  5).  A  collection  of  these  "instantaneous,,  measurements  is  used  to 
statistically  determine  the  parameters  of  interest. 

"instantaneous'1  indicates  a  short  measurement  time  with  respect  to 
the  flow  velocity  variation  of  interest  (i.  e. ,  the  turbulent  frequency).  A 
typical  LV,  measuring  a  100  ft/sec  flow  velocity,  would  have  a  signal 
frequency  of  about  2  MHz,  the  period  of  which  is  0.  5  Msec.  Since  eight 
cycles  are  used  to  measure  the  velocity,  the  measurement  time  is 

4  Msec.  Assuming  an  order  of  magnitude  difference  between  measure¬ 
ment  time  and  the  period  of  the  turbulent  frequency,  the  measurement 
could  be  considered  instantaneous  with  respect  to  a  turbulent  frequency 
of  one  cycle  in  40  Msec,  or  25  KHz. 

Investigations  involving  shear  stress  measurements  (Ref.  6)  indi¬ 
cated  that  the  major  contributions  corresponded  to  frequencies  below 

5  KHz;  therefore,  the  LV  measurements  are  instantaneous  in  many 
flows  of  interest. 
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Figure  3.  Typical  dual  scatter  system. 


3.0  STATISTICS 


The  values  of  the  flow  parameters  determined  with  the  LV  are  the 
results  of  statistical  analysis  of  a  collection  of  individual  measurements. 
Immediately  in  question  is  the  number  of  data  points  necessary  to  ensure 
a  given  accuracy  of  the  parameter  of  interest.  Typically,  estimates  of  the 
required  number  of  data  points  have  been  made  from  approximation  for¬ 
mulas  based  on  the  assumption  of  a  normal  distribution  (Ref.  7).  An 
analysis  not  predicated  upon  this  assumption  was  made,  and  the  results 
were  then  compared  with  values  obtained  from  actual  experimental  LV 
data. 


Consider  a  one -dimensional  flow  having  an  infinite  population  of 
velocities  distributed  about  the  mean  velocity,  p,  with  standard  devia¬ 
tion  <r,  from  which  the  turbulent  intensity  (TI)  =  a/p  can  be  defined.  The 
accuracy  to  which  the  true  values,  p  and  TI,  are  approximated  by  the 
values  v  and  s/v  from  a  smaller  sample  is  desired.  Stated  another  way, 
the  question  is  asked,  what  are  the  variance  of  the  mean  value  (v)  and  the 
variance  of  the  standard  deviation  (s)  of  many  smaller  samples  of  size  N? 
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We  are  interested  in  specifying  a  maximum  tolerance  on  the  mean  ve¬ 
locity  measurement  and  on  the  turbulent  intensity.  Equations  for  the 
statistical  analysis  were  taken  from  Ref.  7,  Chapter  4.  Let  the 
sample  mean  be  indicated  by  v.  The  variance  of  the  mean  of  a  sample 
of  N  measurements  (a£)  is  related  to  the  population  variance  (o^)  by 


o 2  =  a2/ N 


(3.1) 


The  relationship  of  the  standard  deviations  is  then 

<t_  =  o/tJn  (3.  2) 

A  tolerance  on  the  mean  value  of  a  single  sample  (v)  can  be  defined 
as  a  number  (Z)  of  standard  deviations  of  the  mean  (a_)  divided  by  the 
true  (population)  mean  value. 


z  o 

Tol.  <  v)  =  — 


(3.3) 


The  Z  defines  the  confidence  level  of  the  measurement.  For  example, 
a  68-percent  probability  exists  that  a  single  measurement  will  fall  with¬ 
in  la  of  the  mean  (Z  =  1).  Z  =  2  and  3  provide  95-percent  probability 
of  less  than  2a,  and  98 -percent  probability  of  less  than  3a  deviation  from 
the  mean,  respectively.  Relating  to  the  population  parameters,  then 


Tol.  ( v) 


(3.4) 


and  since  the  turbulent  intensity  (TI)  is  the  standard  deviation  of  the 
population  divided  by  the  mean,  then 


Tol.  ( v) 


Z(TI) 

Vn 


It  is  seen  that  the  tolerance  on  the  sample  mean  velocity  is  dependent  on 
the  population  turbulent  intensity. 

For  the  determination  of  tolerance  of  the  turbulent  intensity  we 
introduce  the  r  el -variance  of  a  parameter  (x),  which  is  defined  as 


V 


2 

X 


_2 


x 


(3.5) 
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where  (the  square  of  the  standard  deviation)  is  the  variance  of  x,  and 
x  is  the  mean  value  or  expected  value. 

Letting  s  represent  the  sample  standard  deviation,  the  rel-variance 
of  the  estimated  variance  for  a  simple  random  sample  is 


For  large  N 


\i2  fvi  /S  —  1 

v  *  ~1T 

P  is  called  the  kurtosis  of  the  distribution  and  is  defined  as 


(3.6) 


(3.7) 


where  M 4  =  e  [(v^  -  P)^J  is  the  expected  value  of  the  fourth  moment 
about  the  mean  and  a2  =  c[(v^  -  p)2]  is  the  expected  value  of  the  second 
moment  of  v.  The  value  of  kurtosis  is  always  three  for  a  Gaussian  dis¬ 
tribution. 


The  rel-variance  of  the  estimated  standard  deviation  is  given  by  the 
appr  oximation 


i/2  cv> 


1/4  V2„ 


(3.8) 


The  equation  will  give  good  results  for  moderately  large  values  of  N 
where 


Using  Eq.  (3.7), 


V  2  ^  <  0.3 

V2  *  ft-  1 

S  4N 


(3.  9) 


The  rel-variance  of  the  ratio  of  two  random  variables,  as  in  the 
case  of  the  measured  turbulent  intensity  (s/v)  is,  neglecting  correlation 
terms. 
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V2 


2 

h[i 

or,  in  terms  of  the  relative  turbulent  intensity, 

v(s/»)  =  W  [^T“  +  ™2] 

The  omission  of  the  correlation  term  in  Eq,  (3.  10)  was  experimentally 
justified  as  seen  below  by  the  good  agreement  with  experiment.  In 
general,  if  the  distribution  is  near  normal,  the  correlation  term  is 
small,  and  in  no  case  would  the  discrepancy  be  more  than  ±  V^fe/v)’ 

We  define  the  tolerance  of  the  turbulent  intensity  as 

T  1  /TD  ^  Sample  Std.  Deviations  of  (TI) 

°  ’  Expected  Value  of  (TI) 


(3.  10) 


(3.  11) 


i/2  in*  i/2 

V(s/v)  =  Vs  + 


(\j  —  l 

4N 


ro  /3  -  I 

4N 


which  is  by  definition  ZV 


(s/v) 


or,  finally. 


Tol.  (TI) 


.  z^TpIITn? 


(3.  12) 


It  is  seen  that  the  accuracy  of  the  measurement  of  both  the  mean  velocity 
and  the  turbulent  intensity  is  dependent  on  the  level  of  turbulence  in  the 
flow.  Figure  4  shows  the  sample  size  required  to  ensure  a  given  toler¬ 
ance  on  mean  velocity  for  five  values  of  TI.  A  similar  plot  for  turbulent 
intensity  is  seen  in  Fig.  5,  assuming  a  normally  distributed  population 
(kurtosis  =3). 

An  experiment  was  set  up  to  ascertain  the  validity  of  the  analytical 
results.  A  single  component,  dual  scatter  LV  was  used  to  measure 
velocities  in  a  low-speed  V/STOL  wind  tunnel.  A  turbulence  generator 
was  placed  in  the  flow  to  give  a  moderate  turbulence  level  which  was 
found  to  be  12  percent  upon  data  reduction.  A  large  number  of  data 
points  (10,000)  which  would  theoretically  assure  (with  95-percent  con¬ 
fidence)  a  2 -percent  accuracy  on  TI  for  a  10 -percent  turbulence  level. 
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Tolerance  on  v,  percent 


Figure  4.  Tolerance  on  mean  velocity  versus  sample  size 
for  2 a  confidence  level. 
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Figure  5.  Tolerance  on  turbulent  intensity  versus 
sample  size  for  2 o  confidence  level. 
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were  taken  at  one  position  in  the  flow.  A  histogram  of  the  collection 
showing  number  of  data  points  per  velocity  increment  of  1  ft/sec  is 
seen  in  Fig.  6.  From  this  population,  100  samples  with  N  number  of 
data  points  were  selected  for  several  values  of  N,  using  simple  ran¬ 
dom  sampling  with  replacement.  The  random  sampling  was  used  to 
negate  any  low  frequency  velocity  variations  which  might  have  been 
caused  by  a  slight  drift  in  tunnel  conditions  during  data  collection. 


Velocity,  ft/sec 

Figure  6.  Velocity  distribution,  N  =  10,000. 

For  each  sample  of  N  points,  the  mean  velocity  and  the  turbulent  inten¬ 
sity  were  calculated.  This  produced  100  values  of  v  and  TI.  The 
standard  deviation  of  these  100  values  is  a  measure  of  the  accuracy  one 
could  expect  by  taking  a  single  sample  of  size  N.  To  parallel  the  analyt¬ 
ical  study  the  tolerance  was  defined  to  be  two  standard  deviations  of  the 
parameter  divided  by  the  mean  value  based  on  all  10,  000  data  points. 
This  gives  95-percent  confidence  in  the  measurement. 

2a~  2  ctti 

Tol.  (  v)  =  — -  Tol.  (TI)  =  — - 
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The  comparison  of  the  experimental  and  analytical  values  was  ex¬ 
tremely  good  for  mean  velocity,  as  seen  in  Fig.  7. 


Figure  7.  Tolerance  on  mean  velocity  experiment  versus  theory 
for  turbulent  intensity  =  12  percent. 
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The  comparison  of  the  turbulent  intensity  data  was  very  poor 
assuming  a  normal  distribution  =  3).  The  kurtosis  was  calculated 
for  the  distribution  and  found  to  be  40.  3.  With  j3  =  40,  the  analytical 
values  and  the  experimental  data  compared  well,  as  seen  in  Figure  8. 
For  0  =  40  the  analytical  approximations  used  are  good  for  N  values 
greater  than  433 . 


Figure  8.  Tolerance  on  measurement  of  turbulent  intensity  versus  sample  size 
for  Tl  =  12  percent  and  2a  confidence  level. 
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The  impact  of  this  is  that  for  this  set  of  data  from  a  typical  dual 
scatter  LV,  an  order  of  magnitude  increase  in  the  number  of  data  points 
is  required  to  assure  the  same  accuracy  on  TI  that  one  would  expect  for 
a  normal  distribution. 

An  important  question  is  whether  the  kurtosis  value  can  be  influenced 
by  the  turbulence  source,  turbulence  level,  and  other  flow  character¬ 
istics.  The  answer  is  yes,  and  it  was  also  seen  to  reach  extreme  values. 
To  verify  this  conclusion,  a  second  collection  of  10,  000  data  points  was 
taken  in  the  free  stream  of  the  V/STOL  wind  tunnel  (Fig.  9).  The  turbulent 
intensity  calculated  from  this  collection  was  found  to  be  3.  1  percent. 

The  kurtosis  value  was  found  to  be  1,  533.  The  analytical  prediction 
using  this  value  is  shown  with  the  values  experimentally  obtained  from 
the  population  in  Fig.  10.  Again,  good  agreement  is  seen  for  large  values 
of  N  where  the  approximation  formulas  are  valid.  Of  course,  the  toler¬ 
ance  values  greater  than  100  percent  are  only  for  those  values  of  turbu¬ 
lent  intensity  greater  than  actual.  The  disturbing  prospect  of  a  widely 
varying  kurtosis  indicates  that  in  an  unknown  turbulent  flow,  an  accuracy 
versus  sample  size  prediction  cannot  be  made  prior  to  data  acquisition. 


Velocity,  ft/sec 

Figure  9.  Free-stream  velocity  distribution,  N  =  10,000. 
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Figure  10.  Sample  size  versus  tolerance  on  turbulent  intensity 
for  2 a  confidence  level,  Tl  =  3  percent. 

Laser  velocimeter  turbulent  intensity  measurements  in  a  fully  de¬ 
veloped  channel  flow  with  relative  turbulent  intensities  between  4  and 
6  percent  are  reported  in  Ref.  8.  It  is  shown  there  that  a  slight  varia¬ 
tion  of  kurtosis  (labeled  flatness  factor)  existed  through  the  flow  pro¬ 
file,  but  it  was  always  near  the  value  of  3,  corresponding  to  a  normal 
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distribution.  It  is  believed  that  the  significant  difference  in  the  two 
experiments  is  that  many  of  the  large  particles  were  filtered  out  of  the 
flow  in  the  Naval  Ordinance  Laboratory  (NOL)  study,  indicating  that 
particle  dynamics  may  be  the  primary  cause  for  the  kurtosis  variation. 
Also,  Chauvenet's  criterion,  which  neglects  all  data  over  2or  from  the 
mean,  was  used  in  the  reduction  of  the  NOL  data,  where  no  data  were 
rejected  in  the  present  case.  Data  rejection  is  a  source  of  debate,  and 
of  course  the  two  factors  are  tied  together  in  that  most  of  the  data  points 
rejected  may  be  from  large  particles  lagging  the  flow.  It  is  fair  to 
assume  then  that  for  some  experimental  configurations  (where  sufficient 
flow  filtering  is  possible)  the  problem  may  be  attacked  at  the  data  acqui¬ 
sition  end,  giving  a  reasonable  number  of  data  points  required.  Further 
study  of  the  distribution  shape  factor  is  in  order. 

It  is  important  to  emphasize  that  the  discussion  is  concerned  with 
probabilities  and  tolerances  and  not  actual  error  magnitudes  occurring 
with  a  given  sample  size.  It  can  only  be  stated  that  95  times  out  of  100 
the  values  obtained  from  a  single  statistical  sample  of  size  N  will  fall 
within  the  2a  limit.  Indeed,  68  percent  of  the  time,  the  value  will  fall 
within  half  that  limit,  and  of  course  the  most  probable  value  is  the  mode, 
or  peak,  value. 


4.0  DATA  ANALYSIS 


The  raw  data  output  of  the  LV  is  in  terms  of  the  period  <t)  of  the 
Doppler  signal.  The  desired  parameters  are  in  terms  of  frequency,  or 
1/t.  Two  methods  of  reducing  LV  data  have  typically  been  used.  One 
method  is  to  find  the  mean  period  of  a  collection  of  data  and  from  that  to 
calculate  the  parameters.  The  other  is  to  convert  each  period  measured 
to  a  frequency  and  then  carry  but  the  statistical  manipulations  on  fre¬ 
quencies.  Basic  considerations  would  say  that  if  one  were  interested  in 
velocity 


and  if  the  measurement  error  (At)  of  the  period  was  small,  then  fre¬ 
quency  averaging  would  be  the  proper  way  for  an  unbiased  estimator. 
However,  consider  the  following.  For  a  counting -type  LV,  the  rate 
at  which  scatter  centers  enter  the  probe  volume  is  dependent  on  the 
particle  density  [n(t)],  the  velocity  [v(t)J,  and  the  effective  sampling 
area.  A;  i.  e. , 

N  =  n  |v|  A  (4.  !) 
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The  symbol  N,  representing  particle  flux,  should  not  be  confused  with 
sample  size  from  the  preceding  section.  The  absolute  value  sign  can 
be  removed  if  the  flow  is  in  one  direction  only.  It  is  obvious  that  in  a 
burst-type  situation  (i.  e. ,  where  the  data  are  not  continuous)  a  collec¬ 
tion  of  data  would  be  weighted  toward  the  higher  velocities  in  a  turbu¬ 
lent  flow.  This  is  a  biasing  effect  on  both  the  mean  velocity  and  the 
turbulent  intensity.  An  analysis  of  this  biasing  is  presented  for  the 
case  of  one -dimensional  turbulent  flow. 

The  LV  is  assumed  to  be  operating  in  the  regime  where  the  indi¬ 
vidual  measured  velocities  are  constant  during  the  probe  volume  tran¬ 
sit  time.  That  is,  the  turbulent  eddies  are  much  larger  than  the  probe 
volume  size.  This  ensures  that  each  particle  transit  produces  only 
one  velocity  measurement.  As  individual  velocity  readings  in  turbulent 
flow  have  little  meaning  taken  singly,  average  quantities  are  therefore 
taken. 

A  further  restriction  must  now  be  imposed.  The  following  analysis 
is  based  on  the  assumption  that  the  count  rate  is  totally  determined  by 
the  particle  arrival  rate  in  the  probe  volume.  This  is  true  only  for  the 
case  where 


where  tr  is  the  system  dead  time.  The  net  result  of  large  dead  times  is 
to  reduce  the  bias  for  certain  variable  estimators.  The  analysis  which 
follows  will  set  limits  on  the  biases  which  can  be  expected  from  an  indi¬ 
vidual  realization  LV  system. 


If  Q  is  some  well  behaved  function  of  velocity,  then  the  measured 
average  of  the  quantity  is  given  by 


Qm 


<N  Q> 
<N> 


(4.2) 


where 


Due  to  the  discrete  nature  of  the  counting  process  this  is  true  only 
for  large  numbers  of  particle  events  and  when  the  sample  rate  is  high 
with  respect  to  the  velocity  fluctuations  of  interest.  For  a  more  de¬ 
tailed  discussion  of  the  validity  of  the  assumptions  inherent  in  Eq.  (4.  2) 
the  reader  is  referred  to  Appendix  A  by  D.  O.  Barnett  and  Ref.  9.  To 
determine  the  number  of  events  required  to  assert  that  the  mean  velocity 
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is  accurate  to  within  a  predetermined  accuracy,  the  analysis  of  Section 
3.  0  is  used.  Frequency  averaging  refers  to  the  arithmetic  averaging  of 
the  periods  first  converted  to  frequency.  The  set  of  frequency  samples, 
{f  . . .  f^},  is  related  to  the  set  of  period  samples,  {T^,  . . .  T^j}  by  the 
equation 


The  measured  mean  velocity  as  computed  by  the  two  methods  will 
now  be  derived.  The  results  will  be  given  in  terms  of  the  actual  flow 
parameters,  mean  velocity  and  turbulent  intensity. 

For  both  cases,  period  and  frequency  averaging,  the  number  den¬ 
sity  of  countable  particles  will  be  assumed  to  consist  of  a  mean  value,  n, 
and  a  time -fluctuating  component,  n  ,  such  that 

<n>  =  <n  +  n'>  =  n 

Similarly  for  the  velocity, 

<v>  =  <  v  +  v  '>  =  v 


4.1  MEASURES  OF  THE  MEAN  VELOCITY 

When  a  variable  appears  subscripted  by  m,  followed  by  T  or  f,  the 
quantity  is  a  measured  estimator  of  that  variable  averaged  by  period  or 
frequency,  respectively.  For  example,  is  a  velocity  estimator  ob¬ 
tained  by  frequency  averaging.  From  the  count  rate  model  previously 
introduced,  the  measured  frequency- averaged  mean  velocity  in  terms 
of  flow  parameters  is 


(4.4) 

(4.5) 


_  _  k<nvf>  _  <nv2> 

m  ^  <nv>  <nv>  (4,6) 

Upon  substitution  for  the  number  density  and  velocity  relationships, 
Eqs.  (4.4)  and  (4.5),  and  keeping  only  second-order  terms,  Eq.  (4.6) 
becomes 


<v'2>  <n  'v  '> 

Tj-  +  2  — 

V4  n  v 


1  +■ 


(4.7) 
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A  plot  of  the  percentage  error  as  a  function  of  the  correlation  is  shown 
in  Fig.  11.  Assuming  that  n'  is  independent  of  v' ,  Eq.  (4.  7)  may  be 
written  as 


0  5  10  15  20  25  30  35  40 


Turbulent  Intensity,  percent 
Figure  11.  Fraction  error  in  measured  velocities. 
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This  assumption  is  not  generally  valid.  For  example,  a  mixing  region 
will  show  a  large  correlation  if  the  sources  have  different  particle  con¬ 
centrations.  The  term  <v'^  >/v^  is  just  the  square  of  the  actual  turbu¬ 
lent  intensity.  This  constitutes  a  biasing  of  the  calculated  mean  toward 
higher  velocities,  which  can  be  significant  for  high  turbulent  intensities. 
The  left-hand  side  of  Eq.  (4.  7),  and  consequently  that  of  Eq.  (4.  8),  are 
the  measureable  quantities.  The  remaining  terms  are  yet-unknown 
flow  parameters. 


The  period-averaged  velocity  in  terms  of  flow  parameters  is 


_  A  K  _  K  <n  v> 
vmT  y  <nvT> 


<nv>  =  v/l  *  <dV>N\ 
<n>  \  n  v  J 


(4.9) 


Again  assuming  the  independence  of  n'  and  v' ,  one  finds  the 
bracketed  term  is  unity.  The  velocity  calculated  from  the  periods  is 
found  to  accurately  measure  the  true  mean  velocity  of  the  flow  inde¬ 
pendent  of  the  turbulent  intensity. 


4.2  MEASURES  OF  THE  TURBULENT  INTENSITY 


The  turbulent  intensity,  as  it  would  be  measured  by  the  LV,  can 
also  be  calculated  in  terms  of  the  flow  parameters.  The  frequency- 
averaged  turbulent  intensity  is  given  by 


<  v'2> 


m  f 


<N(y—  Vm{)2> 
<*>  *»f 


<  N  v2> 


(4. 10) 1 


The  measured  velocity  must  be  used  since  the  true  velocity  is  unknown. 

Upon  substituting  for  n  and  v  as  before  and  neglecting  all  terms 
greater  than  second  order,  one  arrives  at 


(TI) 


2  M 

mf 


1  + 


Lv^>/  _  <vf>  _  <nV>\  _ 

t2  \  v2  u; 


<  v'2> 
-2 


+  4 


<n  V>\ 

n  v  / 


<  n  \  '> 
n  v 


(4.11) 
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The  previously  calculated  value  for  the  measured  frequency- 

averaged  velocity,  Eq.  (4.7),  was  introduced  into  Eq.  (4.10).  A  plot  of 

the  error  induced  in  the  frequency -averaged  turbulent  intensity  is  shown 

in  Fig.  12.  .  ,  *  , — 

<  n  v  >  /  n  v 


Figure  12.  Fraction  error  in  measured  frequency-averaged  turbulent  intensity 
with  <n'v'>/  n  v  as  a  parameter. 


The  mean  square  of  the  period  fluctuations  is  now  calculated.  Strictly 
speaking,  this  is  not  a  velocity  turbulent  intensity.  However,  it  does  pro¬ 
vide  an  estimator  for  the  dimensionless  turbulent  intensity,  which  can  be 
obtained  both  theoretically  and  experimentally. 


^2  <T  >mT  <nt2> 


">T 


‘mT 


;  np  2 


N  T 


lil  1-2 


-  1 


mT 


N  T 


mT 


<nV>\ 
n  v  / 


-  1 


(4.  12) 
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Expanding  1/(1  +  — )  into  a  series,  substituting  for  vm<p  and  keeping 
only  second-order  terms,  then 


The  fraction  error  as  a  function  of  turbulent  intensity  is  plotted  in 
Fig.  13.  Assuming  the  correlation  between  particle  and  velocity  fluc¬ 
tuations  is  negligible,  then 


<v'2> 


Again,  as  with  the  mean  velocity,  period  averaging  leads  to  an  unbiased 
estimate. 

<  n'  v' >  /  nv 


Turbulent  Intensity,  percent 

Figure  13.  Fraction  error  in  measured  period-averaged  turbulent  intensity 
with  <nV>/  n  v  as  a  parameter. 

The  equations  for  the  mean  velocity  and  turbulent  intensity  estimators 
can  be  used  to  obtain  the  true  mean  velocity,  turbulent  intensity,  and 
velocity -density  correlation.  Starting  with  Eqs.  (4.  11)  and  (4.  13)  and 
solving  for  <n'v'>/nv  leads  to  a  quartic  equation.  Only  one  of  the  roots 
leads  to  a  physically  meaningful  solution.  The  corresponding  turbulent 


27 


AEDC-TR-74-56 


intensity  and  the  velocity  density  correlation,  when  substituted  into  the 
equations  for  the  mean  velocity,  Eqs.  (4.  7)  and  (4.  9),  should  yield  the 
same  true  velocity. 

Figure  14  shows  the  true  turbulent  intensity  as  a  function  of  the 
period-  and  frequency -averaged  turbulent  intensity.  From  an  experi¬ 
mentally  determined  pair  of  turbulent  intensity  estimators  the  true 
turbulent  intensity  can  be  determined.  Similarly,  the  velocity-density 
correlation  has  been  determined  and  plotted  in  Fig.  15.  Some  care 
should  be  taken  when  the  theory  is  applied  to  specific  problems,  as  the 
assumptions  are  somewhat  restrictive.  Data  should  also  be  free  of  any 
other  broadening  or  biasing  effects  which  may  be  encountered  in  LV 
systems,  such  as  noise-induced  counts.  Again,  the  theory  outlined  is 
valid  only  for  data  rates  which  are  controlled  by  the  particle  arrival 
rates  rather  than  by  the  processor  dead  time,  and  are  high  with  respect 
to  the  velocity  fluctuations. 


Figure  14.  True  turbulent  intensity  as  a  function  of  period-  and 
frequency-averaged  turbulent  intensity. 
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Figure  15.  <nV>/  n  v  as  a  function  of  period-  and 
frequency-averaged  turbulent  intensity. 

5.0  EXPERIMENTAL  DATA 

Laser  velocimeter  data  were  taken  in  a  fre e-jet -mixing  facility  to 
obtain  mean  velocity,  relative  turbulent  intensity,  and  mean  now  angu¬ 
larity.  Hot-wire  anemometer  measurements  were  scheduled  under 
another  project  for  correlation  purposes,  but  that  phase  of  the  project 
was  not  accomplished.  Although  there  is  no  basis  for  comparison  of 
the  LV  data,  it  is  deemed  useful  in  demonstrating  the  technique  and  will 
be  briefly  presented. 

^Facility  operation,  along  with  pitot  probe  and  hot  wire  measure¬ 
ments,  was  accomplished  by  personnel  of  the  Engine  Test  Facility 
Fluid  Dynamics  and  Propulsion  Group  in  Research  Test  Cell  R1A1. 
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The  experiment  comprised  two  concentric  nozzles  arranged  as 
shown  in  Fig.  16.  The  plenum  pressures  were  set  to  give  a  nominal 
2:1  inner  core  to  outer  core  velocity  ratio,  and  the  mixing  of  the  jet 
exhausts  formed  the  turbulent  field. 


Figure  16.  Free-jet-mixing  configuration. 

The  LV  was  a  dual  scatter  system  operated  in  the  forward  scatter 
collection  mode.  It  was  set  up  to  measure  either  of  two  velocity  com¬ 
ponents,  one  at  33.  7  deg.  above  horizontal  and  the  other  at  42.  4  deg. 
below.  The  components  were  not  measured  simultaneously  because  of 
limited  equipment,  but  were  taken  alternately  at  each  location.  This 
should  have  no  effect  on  the  data  since  the  mean  velocity  and  turbulent 
intensity  level  must  be  considered  constant.  At  a  rate  of  six  data  points 
per  second,  a  16 -sec  period  (long  with  respect  to  typical  flow  variations) 
is  required  to  acquire  100  data  points.  This  count  rate  is  much  less 
than  the  rate  of  flow -field  fluctuations. 

The  system  was  mounted  on  a  traverse  which  scanned  radially 
downward  from  centerline.  Radial  profiles  were  taken  at  five  axial 
stations. 

One  hundred  data  bursts  were  processed  for  each  component  at  each 
location.  Under  the  assumptions  of  a  normal  distribution  and  a  10-percent 
turbulent  intensity  this  would  have  assured  (with  95-percent  confidence) 
a  sample  mean  velocity  within  at  least  2  percent  of  true,  and  the  turbu¬ 
lent  intensity  would  not  be  off  more  than  16  percent  (see  Figs.  4  and  5). 
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The  data  were  first  reduced  to  show  mean  velocity  and  standard 
deviation  for  each  velocity  component. 


N 

s 

i=l 


i=l  ‘ 


a  = 


N  -  1 


The  two  velocity  components  were  then  resolved  into  the  mean  velocity 
vector  (V),  the  magnitude  of  which  is 


_ v_\ _ 

(cos  co  s  a  +  sin  Oj  sin  a) 


where 


Tvg  COS  0]  -  Vj  cos  Q£l 

a  =  tan'1  - 2 - — a~ 

vj  sin  02  -  v2  sin 

(The  angles  Q\  and  02  are  the  angles  of  the  measured  components  above 
and  below  earth  horizontal,  respectively. ) 

The  relative  turbulent  intensity  was  then  calculated  in  the  direction 
of  the  measured  components: 

TI  =  vm  and  a2/\  V| 

These  data  are  presented  for  three  axial  stations  in  Fig.  17.  The  turbu¬ 
lence  appears  to  be  nearly  isotropic  in  the  laminar  region  of  the  center^ 
flow,  but  definitely  not  in  the  mixing  region,  where  o^/M  is  twice  ci/|V| 
at  one  point.  With  no  supporting  data  showing  the  shape  of  the  turbulence 
envelope,  no  attempt  was  made  to  adjust  the  data  to  the  flow  coordinate 
system. 

The  data  appear  to  be  internally  consistent.  The  spreading  of  the 
mixing  region,  as  well  as  the  decrease  in  peak  intensity,  is  seen  at 
successive  downstream  stations.  It  is  seen  that  the  peak  occurs  near 
the  middle  of  the  mixing  region.  The  extreme  values,  between  2  and 
20  percent,  appear  reasonable.  Finally,  the  data,  especially  near  center- 
line  at  X  =  3  in.  (Fig.  17a),  show  a  smooth,  fairly  constant  curve  in  that 
region. 
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Radial  Distance,  In. 

a.  Axial  distance  =  3  in.  from  nozzle  exit 


b.  Axial  distance  =  6  in.  from  nozzle  exit 


c.  Axial  distance  =  9  in.  from  nozzle  exit 
Figure  17.  Relative  turbulent  intensity  profiles. 
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Hot-wire  measurements  of  TI  made  for  probe  shakedown  purposes 
were  accomplished  in  nominally  the  same  flow  configuration  (2:1  ratio) 
as  the  LV  measurements.  They  were  limited  to  two  locations  and  were 
oriented  in  a  different  coordinate  system;  however,  their  values  are  of 
interest  for  general  comparison  with  the  LV  data  shown  in  Fig.  17.  On 
centerline  near  the  nozzle  exit  a  value  of  3  percent  for  the  turbulent  in¬ 
tensity  was  obtained  with  the  hot  wire  in  the  200  ft/sec  flow.  This 
agrees  well  with  data  on  one  LV  component  taken  on  centerline  at  each 
of  the  three  axial  stations  shown.  The  other  value,  13  percent,  was 
obtained  in  the  center  of  the  mixing  region.  These  data  fall  between 
those  of  the  two  LV  components  in  all  cases  in  that  region. 


The  total  relative  turbulent  intensity  is  given  as 


(Rel.  TI) 


3 


where  the  subscripts  indicate  the  three  orthogonal  directions  in  the  flow 
coordinate  system.  For  the  purpose  of  calculating  this  property,  it  was 
assumed  that  the  unmeasured  third  component  (03)  was  equal  to  the 
smaller  of  the  other  two.  In  this  case  there  may  be  a  basis  for  assum¬ 
ing  that  it  is  even  much  smaller  since  the  azimuthal  velocity  should  be 
very  small  and  the  energy  gradient  is  radial.  The  absolute  lower  limit, 
where  03  =  0,  and  an  upper  reference,  where  03  =  (a^  +  ct2)/2,  are  shown 
with  these  data  in  Fig.  18  along  with  the  mean  velocity  profile  and  the 
mean  flow  angularity.  Again  the  3  percent  on  the  centerline  and  13  per¬ 
cent  in  the  mixing  region  from  the  hot  wire  compare  well. 


The  mean  flow  angularity  is  referenced  to  earth  horizontal.  It 
appears  that  the  inner  core  jet  was  about  0.  5  deg  above  horizontal.  An 
outward  flow  is  indicated  as  the  angle  becomes  more  negative,  corre¬ 
sponding  to  the  radius  at  which  the  velocity  begins  to  decrease,  i.  e. ,  at 
the  edge  of  the  mixing  region.  In  the  outer  half  of  the  mixing  region  the 
mean  flow  angle  is  radially  inward.  Near -horizontal  flow  is  again  seen 
in  the  low  velocity  outer  jet. 


Pitot  probe  measurements  of  velocity  were  made  concurrently  with 
the  LV  measurements.  The  LV  measurements  were  made  upstream  of 
the  pitot  probe  to  prevent  interference.  The  two  systems  were  mounted 
on  separate  traverse  systems,  and  it  was  evident  that  the  flow  center- 
line  determination  was  different  for  the  two.  Although  the  centerline 
determination  for  the  pitot  data  was  more  accurate,  the  pitot  data  were 
arbitrarily  made  to  coincide  in  radius  with  the  LV  data  for  this  presen¬ 
tation. 
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Radial  Distance,  in. 

Figure  18.  Velocity  and  turbulent  intensity  profiles 
from  LOV  data  showing  flow  angularity. 

A  moderate  discrepancy  (from  4  to  5  percent)  between  the  data  of 
the  two  systems  is  evident.  The  discrepancy  is  believed  to  be  caused 
by  the  data  reduction  technique  used  on  the  LV  data.  Unfortunately,  the 
proposed  new  technique  requires  more  data  points  than  were  taken;  thus 
post-facto  correction  is  impossible.  Two  points  on  the  profile  lent 
themselves  to  the  new  technique.  They  are  shown  in  blackened  symbols  in 
in  the  figure,  giving  0.  3 -percent  agreement  in  the  center  flow  and  about 
2 -percent  agreement  in  the  low  velocity  region. 


6.0  INSTRUMENTAL  DATA  BROADENING 


The  measured  distribution  is  a  convolution  of  the  actual  distribution 
of  velocities  in  the  flow  and  that  of  any  broadening  effects  of  the  system. 
Two  known  systematic  effects  were  investigated  and  found  to  be  negligible. 

There  is  a  slight  variation  in  fringe  spacing  on  the  outer  extremes 
of  the  beam  crossover  region  in  a  dual  scatter  system.  This  is  due  to 
the  Gaussian  intensity  cross  section  of  the  interfering  beams.  For 
typical  beam  intersection  angles  it  is  on  the  order  of  0. 1  percent  and 
is  eliminated  completely  with  an  off-axis  collection  system  since  the 
questionable  region  is  not  seen  by  the  collector  lens  system.  (See 
typical  probe  volume  in  Refs.  10  and  11. ) 
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It  was  also  suspected  that  a  distribution  of  measurements  taken  in 
a  low  signal -to -noise  ratio  (S/N)  environment  would  be  broadened  by  the 
noise  pulses  riding  on  the  signal.  This  possibility  was  investigated 
experimentally  using  a  white  noise  generator  (a  high  gain  amplifier), 
along  with  the  LV  Frequency  Burst  Signal  Synthesizer,  and  the  LV  Data 
Processor,  both  described  in  Ref.  5.  The  noise  signal  was  superim¬ 
posed  on  the  synthesized  LV  signal  before  it  entered  the  data  processor. 
A  signal-to-noise  rms  power  ratio  was  used.  The  rms  signal  amplitude 
was  measured  on  an  oscilloscope,  and  the  noise  amplitude  was  mea¬ 
sured  with  a  broadband  rms  voltmeter  to  set  a  given  S/N.  A  constant 
signal  frequency  of  100  KHz  was  used  to  give  a  basic  period  of  10  psec. 
A  collection  of  at  least  100  data  points  was  recorded  from  the  processor 
for  each  condition.  The  mean  period  (t)  and  the  relative  deviation 
(ct/t)  were  calculated.  The  data,  presented  in  Table  1,  show  a  depend¬ 
ence  on  both  S/N  and  signal  amplitude  for  the  relative  deviation.  In 
the  absolute  extreme  of  poor  S/N,  0.  024,  the  relative  deviation  is  on 
the  order  of  0.  5  percent. 


The  variance  of  a  measured  velocity  distribution  <*m  is  related  to 
the  actual  velocity  distribution  (cr^)  by  cr^  =  +  cr^  where  cr£  is  the 

variance  of  the  noise  distribution;  therefore, 


a 


Table  1.  Data  Broadening  Due  to  Noise 


S/N 

(Power  Ratio) 

Signal 

Amplitude, 

MV 

Noise 

Amplitude, 

MV 

Period - 
Averaged  f , 
jusec 

ot/t, 

percent 

1:1 

7.0 

7.0 

10. 045 

0.  09 

1:1 

3.  5 

3.5 

10. 078 

0.  11 

1:1 

1.  75 

1.  75 

10. 143 

0.  39 

0.2 

7.0 

15.0 

10.043 

0.  12 

0.2 

3.  5 

7.  5 

10.082 

0.  14 

0.  2 

1.  75 

3.  75 

10. 164 

0.  33 

0.  024 

7.0 

42.0 

10.074 

0.48 

0.024 

3.5 

21.  0 

10. 077 

0.  25 

0.024 

1.  75 

11.0 

10. 151 

0.  56 
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Assume  the  extremely  poor  conditions  of  S/N  =  0.  2  and  signal  amplitude  = 
1.  75  mv,  which  would  give  on  =  0.  33  percent.  Assume  also  a  reason¬ 
able  minimum  measured  deviation,  crm  =  3.0  percent.  For  this  "worst 
case"  condition,  a  would  be  2.  98  percent  for  a  0.  67 -percent  error  in 
the  measured  value.  A  well  designed  LV  should  exceed  a  S/N  of  unity 
and  a  7-MV  signal  amplitude,  whereby  the  measured  deviation  error 
should  be  less  than  0.  06  percent.  This  is  not  a  significant  contribution 
to  the  turbulent  intensity  measurement. 


7.0  SUMMARY  AND  CONCLUSIONS 


The  LV  system  is  well  suited  for  direct  measurements  of  turbulent 
intensity.  All  of  the  systematic  broadening  effects  are  insignificant  with 
respect  to  typical  flow  turbulent  intensities  of  3  percent  and  above. 

It  was  found  that  a  much  larger  quantity  of  data  than  previously  ex¬ 
pected  is  required  to  specify  a  given  accuracy  on  a  measurement  of  TI. 
This  is  true  because  the  typical  data  obtained  were  not  normally  distri¬ 
buted  and  exhibited  large  values  for  the  kurtosis.  The  analytical  re¬ 
lationship  defining  required  sample  size  was  verified  for  large  N  values 
by  two  sets  of  experimental  data.  It  has  been  shown  that  the  kurtosis  is 
affected  slightly  by  the  turbulence  source.  If  it  can  feasibly  vary  widely, 
then  an  estimate  of  kurtosis  will  be  needed  before  an  estimate  of  sample 
size  can  be  made  for  each  location  in  a  turbulent  field.  The  other  alterna¬ 
tive  would  be  to  use  a  pessimistic  value  of  kurtosis  and  accept  the  result¬ 
ing  tolerance  on  the  measurements. 

Indications  are  that  kurtosis  is  more  affected  by  other  system  or 
flow  characteristics  such  as  large  particle  dynamics,  which  could  cause 
a  variation  in  the  low  velocity  tail  of  the  distribution.  If  this  is  the  case, 
it  may  be  possible  to  discriminate  against  such  data  points  by  flow  filter¬ 
ing  or  electronic  means,  and  to  decrease  the  required  sample  size. 
Further  studies  are  needed  in  this  area  defining  the  significant  effects 
on  the  characteristics  of  the  measured  distribution. 

It  was  found  that  both  frequency _aver aging  and  period  averaging  pre¬ 
sented  a  biased  estimator  of  TI  and  V.  The  system  of  equations  de¬ 
veloped  to  predict  the  biases  forms  a  basis  for  generalization  to  more 
complex  geometries.  The  period -averaged  mean  velocity,  in  general, 
is  least  biased  and  is  independent  of  turbulent  intensity,  when  sample 
rates  are  very  high.  However,  for  most  practical  applications  this  bias 
is  small,  as  sample  rates  are  controlled  by  instrument  deadtime. 
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APPENDIX  A 

BIAS  EFFECTS  OF  INDIVIDUAL  REALIZATION  VELOCIMETERS, 

by  D.  O.  Barnett 


The  turbulence  at  a  point  in  an  incompressible,  isothermal  flow 
can  be  totally  characterized  if  the  velocity  history,  V(t),  is  known.  In 
LV  practice,  however,  only  certain  velocity  moments  are  obtained 
which,  together  with  basic  assumptions  concerning  the  velocity  distri¬ 
butions,  adequately  describe  the  process  of  interest. 

If  Q  =  Q(V),  then  the  quantity  one  wishes  to  approximate  with  an 
LV  system  is  the  temporal  mean  of  Q,  which  may  be  defined  as  follows: 

<Q>  =  I  ;T  Q(V)dt  (A-i) 

1  o 

Data  taken  with  an  individual  realization  LV  are  not  continuous, 
and  consequently  Eq.  (A-l)  must  be  approximated  as 

(A-2) 

where  {Q}  is  now  an  estimator  for  <Q>,  the  desired  temporal  average  of 

Q. 


While  the  individual  values  are  accurately  known,  the  correspond¬ 
ing  sample  intervals,  which  minimize  bias,  require  a  general  model  of 
the  sampling  process.  Obviously  a  deterministic  prediction  of  the  exact 
arrival  times  at  the  probe  volume  cannot  be  made.  However,  if  in  the 
vicinity  of  the  probe  volume  the  number  density  is  constant,  then  the 
time  interval  between  successive  samples  for  one -dimensional  flow  can 
be  written  as 


A  1  1 

Ati  =  <v.>  *  nA 


(A-3) 


where  <V^>  is  the  time-averaged  fluid  velocity  over  the  interval,  n  is 
the  particle  number  density,  and  A  is  the  probe  area..  Correlations 
between  the  velocity  and  particle  density  have  not  been  included,  but  it 
should  be  noted  that  Eq.  (A-3)  implies  that  there  is  a  strong  correlation 
between  velocity  changes  in  the  probe  volume  and  the  velocity  immediately 
upstream  of  the  probe  volume. 
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Upon  introduction  of  Eq.  (A-3)  into  Eq.  (A-2)  and  noting  that 


T  =  X  At; 
i=i  1 


(A-4) 


one  solves  for  the  temporal  mean  of  Q  thus: 


X  0;/<  v>i 

)Q  i  =  — - 

v  M 

X  1/<V>. 
i=i  1 


(A-5) 


The  conditions  under  which  Eq.  (A-5)  yields  an  unbiased  estimate 
of  <Q>  will  now  be  determined. 

Expanding  V(t)  in  a  series  of  the  form 


V(t)  =  <V>  +  ^X  VN  cos  («Nt  -  Bn) 


(A-6) 


yields  an  expression  which  is  valid  over  the  entire  sampling  time,  T. 

By  integration  (Eq.  A-6)  over  a  sample  interval  Atj,  the  mean  value  can 
be  obtained. 


<v>i 


<v>  + 


X 

N=1 


VN 


sin  (<y 


(A-7) 


Case  I:  Sample  Frequency  Less  Than  Minimum  Flow  Frequency. 


When  UNAti  »  1  (that  is,  when  the  sampling  frequency  is  much  less 
than  the  lowest  frequency  of  the  flow  oscillations  of  interest),  <V>j  -*<V>. 
For  this  case  the  time  interval  of  the  sample  is  essentially  constant,  and 
the  temporal  of  Eq.  (A-5)  and  the  arithmetic  average  are  identical.  The 
arithmetic  average  is  defined  as 


M 


Q  =  1/M  X  Q; 
i=  1 


(A-8) 


It  follows  that  the  mean  velocity  over  the  entire  time  interval,  T,  also 
is  given  by  the  arithmetic  mean  of  the  sample  velocities,  Vit  While 
this  procedure  does  not  prove  the  equivalence  of  the  arithmetic  and 
temporal  averages,  it  should  be  noted  that  for  sample  frequency  low 
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with  respect  to  the  minimum  significant  frequency  component,  any 
correlation  between  the  instantaneous  velocity  and  the  sampling  rate 
tends  to  be  destroyed.  Since  it  is  equally  probable  that  during  any 
sample  interval  the  sampled  velocity  will  be  greater  or  less  than  the 
mean  value,  the  arithmetic  mean  should  provide  an  unbiased  estimator 
of  the  mean  velocity. 


Case  II:  Sample  Frequency  Greater  Than  the  Maximum  Flow  Frequency 

If  the  sample  frequency  is  large  compared  to  the  frequency  of  the 
flow  oscillations,  then  W]\jAt-[  «  1. 


For  this  case,  the  individual  sample  accurately  represents  the 
velocity  during  the  time  interval.  In  essence,  the  flow  appears  station¬ 
ary  during  the  sample  interval,  and  <V>^  -  Vj.  Furthermore,  for  large 
sample  sizes  the  finite  sum,  Eq.  (A-2),  becomes  an  excellent  approxi¬ 
mation  to  the  temporal  mean,  Eq.  (A-l).  Equation  (A-2)  reduces  to  the 
harmonic  average  of  the  individual  realizations  of  the  velocity  thus: 


{VI 


i  M 

M  1  VVi 
M  ;=i  1 


K 


(A-9) 


where  K  is  a  constant  (a  function  of  the  LV  geometry  and  wavelength),  and 
the  period  is  given  by  =  1/f^. 

It  is  thus  Eq.  (A-8)  which  forms  the  basis  for  the  statistical  bias  of 
the  LV  data.  The  interpretation  of  the  harmonic- averaged  velocity  data 
as  a  "period-averaged"  velocity  is  evident. 


One  can  show  that  the  arithmetic  average  of  the  frequency  yields 
higher  mean  values  than  the  harmonic -averaged  result.  This  may  be 
demonstrated  by  letting 


fi  =  f  (A-10) 

Upon  substituting  into  the  expression  for  the  harmonic  mean,  Eq.  (A-8), 
and  using  the  expansion  for  (1  +  x)"1, 

V  =  U  =  IV!  1  (-l)rF/7r  /A  „  % 

r=o  (A- 11) 

for  f '  <  f.  For  symmetric  velocity  distributions  and,  hence,  frequency 

distributions,  the  odd  order  moments  (i.  e. ,  r  =  1,  3,  5  .  .  . )  yield 
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vanishing  central  moments,  establishing  the  inequality  of  the  arithmetic 
and  temporal  averages. 

The  analysis  indicates  that  for  one -dimensional  flow  through  the 
probe  volume  of  an  LV  system,  it  is  possible  to  obtain  arithmetic 
averages  which  are  biased  with  respect  to  the  temporal  mean  velocity 
of  the  fluid.  If  the  sampling  frequency  is  much  higher  than  the  frequency 
of  the  velocity  fluctuations,  the  bias  can  be  as  large  as  predicted  by  the 
analysis.  On  the  other  hand,  sampling  at  low  frequency  will  eliminate 
the  fundamental  bias.  The  situation,  however,  will  be  mixed  in  turbu¬ 
lent  flow  where  many  harmonics  will  occur. 

For  those  sample  rates  where  bias  is  predicted,  the  bias  in  the 
mean  velocity  can  be  eliminated  by  harmonic  averaging  of  the  individual 
realizations  of  velocity. 
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